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I .  INTRODUCTION 


There  has  been  a  long-standing  requirement  for  a  detailed  description 
of  the  mass,  velocity  and  size  distributions  of  fragments  created  during 
projectile  impact  on  armor  targets.  This  information  is  required  for  an 
interpretation  of  the  damage  a  fragment  cloud  could  cause  a  secondary 
target. 

Previously,  information  about  fragment  debris  has  only  been  obtain¬ 
able  through  expensive  experimentation,  but  within  the  last  two  decades, 
there  has  been  substantial  progress  in  the  development  of  an  active  frac¬ 
ture  model  (NAG/FRAG) 1-3  that  is  essentially  sufficient  for  the  prediction 
of  exactly  the  fragment  information  desired. 

The  complete  fracture  model  is  based  on  the  nucleation  and  growth  of 
voids  and  subsequent  void  coalescence,  i.e.,  features  of  failure  that  are 
known  to  occur  with  the  continued  application  of  intense  loads.  In  keep¬ 
ing  with  experimental  observations,4  fragmentation  is  treated  as  the 
coalescence  of  voids  and  is,  therefore,  dependent  on  nucleation. 

Through  investigations  of  several  materials  including  iron,2  copper,5 
beryllium6  and  apolycarbonate,  Lexan, 7  it  has  been  discovered  that  the 
nucleation  rate  N  of  voids  in  many  materials  is  exponentially  dependent 
on  the  applied  stress  a  and  satisfies  the  relation 


Moss  and  C.  M.  Glass 3  " Some  Microscopic  Observations  of  Crooks  Devel¬ 
oped  in  Metals  by  Very  Intense  Stress  Waves, "  Ballistio  Research  Labora¬ 
tories  Technical  Note. No.  1312,  April  1960.  (AD  #237943) 

2T.  Barbee,  L.  Seaman ,  and  R.  C.  Crewdson ,  " Dynamic  Fracture  Criteria  of 
Homogeneous  Materials, "  Air  Force  Weapons  Laboratory  Technical  Report 
No.  AFWL-TR- 70-99,  November  1970. 

3 L.  Seaman,  D.  R.  Curran ,  and  D.  A.  Shookey,  " Computational  Models  for 
Ductile  and  Brittle  Fracture, "  J.  Appl.  Phi) 8. ,  Vol.  47,  No.  11, 

November  1976,  pp.  4814-4826. 

4D.  A.  Shockey,  L.  Seaman,  D.  R.  Curran,  P.  S.  DeCarli,  M.  Austin ,  and 
J.  P.  Wilhelm,  "A  Computational  Model  for  Fragmentation  of  Armor  Under 
Ballistic  Impact, "  Ballistic  Research  Laboratories  Contract  Report 
No.  222,  April  1975.  (AD  #B004672L) 

5L.  Seaman,  T.  W.  Barbee,  Jr.,  and  D.  R.  Curran,  " Dynamic  Fracture 
Criteria  of  Homogeneous  Materials,  "  Air  Force  Weapons  laboratory 
Technical  Report  No.  AFWL-TR-7 1-156,  Deoember  1971. 

eD.  A  Shockey ,  L.  Seaman ,  and  D.  R.  Curran,  "Dynamic  Fracture  of 
Beryllium  Under  Plate  Impaot  and  Correlation  with  Electron  Beam  and 
Underground  Test  Results ,"  Air  Force  Weapons  Laboratory  Technical 
Report  No.  AFWL-TR-7 3-12,  June  1973. 

7D.  R.  Curran,  D.  A.  Shockey ,  and  L.  Seaman,  "Dynamic  Fracture  Criteria 
for  a  Polycarbonate,"  J.  Appl,  Phys.  Vol .  44,  No.  9,  September  1973, 
p.  4025. 
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N  ■  e 


(1) 


where  o  „  is  the  threshold  stress  for  the  nucleetion  of  voids,  and  ft  and 
no  o 

Oj  are  experimentally  determined  material  parameters. 


In  addition,  the  fracture  model  includes  equations  for  the  growth 
and  coalescence  of  voids,  and  these  equations,  in  turn,  depend  on 
independent  sets  of  material  parameters. 


All  the  material  parameters,  including  onQ,  are  usually  established 

simultaneously  with  a  trial-and-error  procedure  -  one  in  which  a  set  of 
parameters  is  assumed  and  used  in  a  computation  to  predict  the  total 
cracking  from  specific  stress  histories.  Then,  the  cracking  of  samples 
actually  subjected  to  these  loading  conditions  is  determined  experi¬ 
mentally  and  compared  with  the  computed  results.  The  procedure  is 
repeated  until  there  is  agreement,  which  implies  the  correct  parameters 
have  been  found,  or  until  continued  disagreement  suggests  the  fracture 
model  is  being  used  with  inappropriate  nucleation  and  growth  functions. 
When  the  approach  is  properly  applied,  sufficient  test  data  is  used  to 
ensure  unique  results  for  the  range  of  conditions  investigated. 


While  all  the  parameters  can  be  found  by  this  trial-and-error  pro¬ 
cedure,  their  determination  is  substantially  simplified  if  onQ  is  found 

by  an  independent  method.  Such  an  independent  evaluation  was  attempted 
and  is  the  subject  of  this  report. 


II.  PROCEDURE 

The  material  used  for  these  tests  was  MIL  SPEC  12560-B  rolled 
homogeneous  steel  armor  (12.7  mm  thick  plate)  with  the  composition  given 
in  Table  I.  Mechanical  properties  of  this  same  plate  of  steel  have  been 
characterized  in  parallel  research  efforts  to  furnish  consistent  data 
for  this  material. 8-10 


8J?.  F.  Benck,  "Quasi-Static  Tensile  Stress  Strain  Curve  i — II,  Rolled 
Homogeneous  Armor,"  Ballistic}  Reeearoh  Laboratories  Memorandum  Report 
2703,  November  1976.  (AD  #B016015L) 

9R.  F.  Benak  and  J.  L.  Robitaille,  " Tensile  Stress-Strain  Curves — III, 
Rolled  Homogeneous  Armor  at  a  Strain  Rate  of  0.42  s~^ , "  Ballistic 
Researoh  Laboratory  Memorandum  Report  2760,  lime  1977.  (AD  #A041560) 

10ff.  E.  Hauver ,  " The  Alpha-Phase  Hugoniot  of  Rolled  Homogeneous  Armor," 
Ballistic  Research  Laboratories  Memorandum  Report  2661,  August  1976. 
(AD  #B012871L) 
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TABLE  I 

COMPOSITION  OF  STEEL,  WEIGHT 

PERCENT 

c 

Mn  P  S 

Si 

Ni 

Cr 

0.22 

0.26  0.001  0.015 

0.19 

3.15 

1.06 

Cu 

V  Mo  A1 

B 

Ti 

0.1 

0.01  0.15/0.30  0.03 

N.D. 

N.D. 

The  microstructure  of  the  material  was  bainitic  with  inclusions  as 
shown  in  Fig.  1.  The  Rockwell  "C"  hardness  was  37,  and  the  material 
will,  therefore,  be  referred  to  in  the  following  as  RHA  (RC-37) . 


Figure  1.  Rolled  homogeneous  steel  armor.  Nital 
etch;  magnification  500. 
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It  was  considered  possible  that  ono  could  be  determined  independently 

of  the  NAG/FRAG  model  if  the  density  of  voids  created  at  several  stresses 
above  the  nucleation  threshold  stress  were  plotted  versus  the  impact 
compressive  stress,  and  the  curve  was  extrapolated  to  the  stress  where  no 
voids  would  form,  i.e.,  to  the  no-damage  condition. 

For  such  a  procedure  to  work,  the  need  to  account  for  special  details 
normally  taken  i..to  account  by  the  NAG/FRAG  analysis  must  be  overcome. 
Specifically,  the  reduction  in  stress  amplitude  by  elastic  wave  reflec¬ 
tions  from  a  plastic  wave  front,  as  occurs  when  plate  impacts  are  used 
to  cause  fracture,  and  the  changes  in  stress  associated  with  unloading  at 
void  interfaces  must  either  be  treated  independently  or  be  reduced  to 
negligible  effects. 

The  approach  used  was  to  try  to  render  the  elastic  wave  reflections 
from  the  plastic  wave  and  the  unloading  at  void  interfaces  negligible  by 
working  within  a  stress  range  just  above  the  nucleation  threshold  stress. 

There  should  be  only  a  weak  wave  reflected  from  an  encounter  of 
elastic  and  plastic  waves  when  the  amplitude  of  the  plastic  wave  is  only 
slightly  greater  than  the  Hugoniot  elastic  limit.11  Furthermore,  the 
average  effect  of  unloading  at  void  interfaces  should  be  trivial  since 
only  a  few  small  voids  would  be  generated  at  low  stresses.  Certainly, 
the  "reflection"  effect  should  go  to  zero  at  the  elastic-plastic  transi¬ 
tion  stress  (EPTS) ,  and  the  "void  interface"  effects  should  go  to  zero 
at  the  nucleation  threshold  stress. 


It  has  been  observed  in  previous  fracture  tests  that  the  nucleation 
threshold  stress  is  often  approximately  the  same  as  the  EPTS,  and  without 
exceeding  it,  there  would  be  little  chance  of  breaking  a  sample.  Hence, 
the  EPTS  was  taken  as  the  lower  bound  on  the  range  of  shock  stresses  of 
interest,  and  the  upper  bound  was  just  one  GPa  higher. 

Examples  of  fracture  from  impact  loading  in  this  stress  range  were 
created  with  parallel  plate  impacts  accomplished  with  a  light-gas  gun. 

The  maximum  compressive  stress  from  such  an  impact  is  given  by10,12 


0=0  +  0 
e  p 


(2) 


where  ofi  is  the  amplitude  of  the  elastic  precursor  that  leads  the  plastic 
part  of  the  stress  wave,  and  o^  is  the  difference  between  the  maximum 


11T.  W.  Wright ,  "Ultrasonic  Probing  of  Plastic  Waves,"  J.  Appl.  Phys., 
Vol.  39 ,  No.  12,  November  1968,  pp.  5740-5745. 

l2W.  H.  Rice ,  R.  G.  McQueen,  and  J.  M.  Walsh,  "Compression  of  Solids  by 
Strong  Shock  Waves, "  Solid  State  Physics,  Vol.  6,  F.  Seitz  and 
D.  Turnbull,  Eds.,  Academic  Press,  New  York,  1958,  pp.  1-63. 
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Figure  2.  (a)  Lagrangian  plot  of  the  location  of  the  elastic  and 

plastic  wave  fronts  versus  time,  (b)  Maximum  tensile 
stress  duration  versus  distance  through  the  thickness 
of  the  plate-shaped  sample. 


It  has  been  found  experimentally  that  the  elastic-plastic  transition 
stress  (EPTS)  of  shock  waves  emanating  from  planar  impact  interfaces 
decreases  in  amplitude  with  travel.  *3  In  the  case  of  the  12.7  mm  thick 
samples  used,  the  shock  wave  traveled  19.1  mm  to  the  zone  of  maximum 
tensile  stress  duration  which  was  roughly  where  maximum  damage  should  be 
expected.  After  a  shock  wave  travels  19.1  mm  in  RHA  (RC-37) ,  the  EPTS 
drops  to  about  1.9  GPa. 10  Hence,  this  was  selected  as  the  lower  bound 
on  the  impact  pressures  of  interest. 

A  summary  of  the  impact  velocities  and  the  resultant  maximum 
compressive  stress  levels  is  given  in  Table  II. 


130.  M.  Butcher  and  D.  E.  Munson,  " The  Application  of  Dislocation  Dynamics 
to  Impact-Induced  Deformation  Under  Uniaxial  Strain ,  "  Dislocation 
Dynamics,  A.  R.  Roeen field ,  G.  T.  Hahn ,  A.  L.  Bement,  Jr.  and  R.  I. 
Jaffee ,  Eds.,  McGraw-Hill  Book  Co.,  Hew  York,  1968,  pp.  591-607. 
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TABLE  II 


TEST  CONDITIONS  USED  TO  ESTABLISH  a 
_  I  no 


IMPACT  VELOCITY 

MAXIMUM  COMPRESSIVE  STRESS 

mm/usec. 

GPa 

0.0954 

2.08 

0.1178 

2.49 

0.1392 

2.88 

In  this  work,  the  impact  velocities  were  measured,  and  the  stresses 
listed  in  Table  II  were  determined  with  Eqs.  2  through  4  and  the  Hugoniot 
relation10 


U,  =  4.51  +  1.43  u 
P  P 


(5) 


III.  RESULTS 

Observation  of  the  recovered  and  sectioned  samples  revealed  they 
were  all  partially  broken  by  the  impacts  and  that  the  nature  of  the 
failure  was  the  same  in  all  the  samples.  However,  the  voids  in  the  low- 
pressure  test  were  confined  entirely  to  a  narrow  region  at  the  midplane  of 
the  plate. 

Microscopic  observation  of  the  sectioned  samples  revealed  the  break¬ 
ing  process  was  one  in  which  (1)  inclusions  cracked  or  became  .separated 
from  the  matrix,  (2)  "hair-line"  cracks  extended  into  the  matrix  from  the 
damaged  inclusions  -  while  under  load,  these  cracks  were  presumably  open, 
(3)  the  "hair-line"  cracks  extended  from  inclusion  to  inclusion  and  (4) 
cracks  on  slightly  separated  planes  coalesced  by  shear  deformation  on 
surfaces  connecting  the  ends  of  the  cracks. 

Clearly,  there  are  several  distinct  stages  in  the  failing  process, 
and  nucleation  can  be  interpreted  in  several  ways.  For  example,  a  stable 
nucleus  could  quite  naturally  be  envisioned  as  a  cracked  inclusion.  Here, 
however,  nucleation  was  associated  with  the  development  of  the  "hair¬ 
line"  cracks  that  extended  from  inclusions  into  the  matrix.  Hence,  cracks 
were  not  intentionally  counted  until  they  had  reached  this  stage  of 
development. 

Identifying  and  counting  cracks  in  this  way  is  a  useful  way  to  keep 
track  of  the  damage  to  RHA  because  the  "hair-line"  cracks  must  nucleate 
and  grow  before  the  material  can  separate  into  pieces. 
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Actually,  the  connection  of  voids  by  plastic  shear  is  the  last 
feature  of  the  breaking  process,  but  the  heating  on  these  shear  zones  by 
plastic  working  results  in  weakened  material  on  which  further  separation 
can  occur  relatively  easily.  It  is  convenient,  therefore,  to  identify 
this  stage  of  failure  with  coalescence  and  to  treat  it  with  a  recently 
developed  fragmentation  analysis.4  Hence,  voids  connected  by  shear  were 
counted  as  two  cracks  when  separated  by  approximately  one-third  the 
distance  between  the  faces  of  the  most  open  crack. 

Eq.  (1),  if  applicable,  implies  that  nucleation  continues  as  long  as 
a  stress  greater  than  o^o  is  applied.  Furthermore,  it  is  clear  from 

Fig.  2  that  the  load  duration  depended  on  the  location  in  the  plate. 
Hence,  only  voids  from  the  central  regions  of  each  plate  were  counted. 

For  the  1.93  GPa  test,  this  was  a  region  0.06b  cm  (0.026  in.)  wide 
centered  about  the  maximum  damage  in  the  plate.  For  the  other  tests, 
this  was  a  region  0.127  cm  (0.05  in.)  wide,  also  centered  about  the 
maximum  damage.  This  ensured  the  voids  from  each  test  were  initiated 
over  equal  time  intervals  and  allowed  the  use  of  data  from  the  low 
pressure  test  -  an  advantage  in  reducing  the  extent  of  an  extrapolation 
to  the  nucleation  threshold  stress. 

A  graph  of  the  total  number  of  voids  counted  in  these  central 
regions  versus  the  maximum  compressive  stresses  imposed  on  the  samples 
is  shown  in  Fig.  3. 


STRESS  (GPo) 


Figure  3.  Stress  dependence  of  void  densities  in  the  regions  of  maxi¬ 
mum  damage  -  approximately  the  midplane  of  each  sample. 
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IV.  DISCUSSION 


It  was  not  known  before  examining  the  results  if  the  stress  range 
selected  for  this  work  was  sufficiently  near  the  nucleation  threshold 
stress  to  make  the  impact  compressive  stresses  a  good  approximation  of 
the  maximum  tensile  stresses  developed.  This  was,  therefore,  checked  by 
computing  the  maximum  tensile  stresses  with  the  one-dimensional  wave- 
propagation  computer  program  PUFF  and  the  brittle  fracture  subroutine 
BFRACT  2.4»14  Such  a  computation  automatically  accounts  for  the  elastic- 
plastic  wave  interactions  as  well  as  the  development  of  voids. 

The  occurrence  and  development  of  voids  is  treated  simultaneously 
with  the  loading  history  which  is  the  only  way  predictions  of  the  process 
can  be  made  because  failure  and  loading  are  interdependent  processes. 

The  failure  (development  of  cracks)  that  occurs  and  the  corresponding 
adjustment  in  the  load  (stresses)  are  based  on  the  volume  change  AV 
imposed  on  the  sample.  It  is  assumed  that  this  change  is  a  composite  of 
the  volume  change  by  strain  in  the  solid  AVg,  volume  change  due  to 

nucleation  AVn  and  volume  change  by  void  growth  and  expansion  AVg.  When 

expressed  in  these  terms,  the  volume  change  is  given  by4 

AV  =  AV  +  AV  +  AV  .  (6) 

s  n  g 

The  change  in  solid  volume  AVg  during  the  computation  time  increment 

is  related  to  the  pressure  in  the  solid  Pg  which  is  obtained  from  the 
Mie-Griineisen  relation4 


P 

s 


=  C 


MVso  + 


V 

so 


rE 

+  AV 

s 


(7) 


where  C  is  the  bulk  modulus,  r  is  the  Griineisen  ratio,  E  is  the  internal 
energy  and  the  subscript  "o"  refers  to  conditions  at  the  beginning  of  the 
time  increment. 


The  increment  of  nucleated  void  volume  AVn  is  determined  from  the 

product  of  the  number  of  voids  nucleated  during  the  time  increment  (Eq.  1) 
and  the  volume  of  the  voids.  The  growth  increment  AV  is  based  on  both 

o 

the  increase  in  crack  radius  and  the  elastic  opening  of  the  crack. 

Together,  the  three  terms  Av  ,  AV  and  AV  are  related  in  complex  ways 

s  n  § 

to  the  applied  stress  so  that  the  stress  cannot  be  obtained  directly 
for  the  imposed  strain  increments. 
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The  solution  for  stress  is  obtained  by  an  iteration  procedure  in  which 
a  stress  tensor  is  first  estimated.  Then,  the  three  increments  on  the 
right  side  of  Eq.  6  are  obtained  directly.  The  sum  of  these  increments 
is  AV  ,  an  approximation  of  AV.  Repeated  stress  estimates  are  made  with 

a  regula  falsi  method  until  AV&  is  sufficiently  close  to  AV.  Results  are 

made  reasonably  independent  of  step  size  by  requiring  that4 

(AV  -  AVayVs  <  2  X  10'5  ,  (8) 

where  AV  is  one  of  the  iterative  estimates  of  the  volume  change.  It  has 

a 

been  found  that  when  this  inequality  is  satisfied,  results  are  independent 
of  step  size  to  a  precision  of  0.003  GPa. 

When  applying  the  BFRACT  mod,el,  it  is  important  to  represent  the 
crack  nucleation  rate  correctly.  This  was  assumed  to  be  given  by  Eq.  1 
which  is  consistent  with  the  material  behavior  if,  for  constant  load 
durations.  In  N  versus  o  is  a  straight  line.  The  crack  densities  shown 
in  Fig.  3  were  obtained  from  regions  of  maximum  damage  for  which  the  load 
duration  was  essentially  the  same  in  each  test,  and  when  the  data  is 
replotted  as  in  Fig.  4,  it  is  clear  that  the  material  behavior  is  repre¬ 
sented  very  well  by  Eq.  1.  Hence,  the  method  of  calculating  the  tensile 
stresses  was  reasonable. 


Figure  4.  Damage  caused  by  tensile  stress  waves  at 
the  midplane  of  samples. 
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The  computed  tensile  stresses  are  shown  in  Fig.  3,  and  it  is  evident 
that  they  are  less  than  the  impact  pressures.  Furthermore,  the  differ¬ 
ence  in  these  stresses  increases  with  the  intensity  of  the  load.  Hence, 
an  extrapolation  of  the  compressive  stresses  to  the  no-damage  condition 
would  give  a  substantially  incorrect  value  for  onQ. 

Presumably,  the  compressive  stresses  could  still  be  used  in  an 
extrapolation  to  the  no-damage  condition  provided  useful  data  were 
obtained  from  tests  conducted  at  stresses  below  2  GPa.  However,  experi¬ 
mental  evidence  suggests  this  would  be  impractical  for  RHA  (RC-37) 
because  the  cracks  formed  are  so  widely  spaced  when  maximum  stresses  are 
in  this  range  that  it  would  be  difficult  to  determine  meaningful  average 
crack  densities. 

Hence,  it  is  concluded  that  it  is  not  possible  to  use  an  extrapola¬ 
tion  independent  of  an  "active"  fracture  analysis  to  determine  o  for 
RHA  (RC-37).  n0 

It  is,  however,  useful  to  extrapolate  the  crack  density  against 
tensile  stress  to  the  no-damage  condition  because  the  nucleation  threshold 
stress  established  with  this  procedure  necessarily  corresponds  to  the 
stress  where  failure  begins.  In  contrast,  the  "trial-and-error"  procedure 
used  to  determine  the  nucleatioi.  and  growth  material  parameters  is  based 
on  the  best  fit  of  observable  fracture  data,  but  these  are  developed  with 
loading  conditions  that  do  not  include  the  stress  where  fracture  begins. 

The  suggested  extrapolation  is  shown  in  Fig.  3,  and  the  nucleation 
threshold  stress  on®  corresponding  to  the  no-damage  condition  is  1.65  GPa. 

tg 

The  result  of  a  trial-and-error  determination,  crno  ,  was  1.8  GPa. 

tg 

The  two  results  are  different  because  in  determing  the  density  of 

cracks  with  radii  less  than  0.15  mm  were  underestimated  because  the  size 

distribution  curve  (In  N  vs.  R)  was  nonlinear  in  the  size  range 

c 

0  <  R  <  0.60  mm  and  more  complicated  than  could  be  described  with  the 
model  at  hand.  Cracks  with  radii  between  0.15  and  0.60  mm  were,  however, 
represented  correctly. 

The  result  of  the  extrapolation  should  only  be  in  error  to  the  extent 
that  cracks  with  radii  in  the  range  0  <  R  <  0.15  mm  were  ignored  while 
computing  the  tensile  stresses  with  the  PUFF-BFRACT  computer  code.  This 
error  should  be  small  because  these  voids  contribute  little  to  the  total 
void  volume.  Hence,  the  actual  nucleation  threshold  stress  of  1.65  GPa 
should  be  reasonably  accurate. 

As  already  indicated,  the  EPTS  at  the  midplane  of  the  samples  was 
1.9  GPa.  Thus,  it  appears  that  cracks  nucleated  at  stresses  below  the 
yield  stress  of  the  material.  This  could  be  due  to  moderate  stress 
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concentrations  at  inclusions  or  the  method  that  has  been  used  to  deter¬ 
mine  the  EPTS.10 

The  EPTS,  as  reported  in  Reference  10,  was  determined  by  extending 
the  straight  line  portions  of  the  stress  wave  profile  (free  surface 
velocity  vs.  time)  on  adjacent  sides  of  the  transition.  The  intersection 
of  these  straight  lines  was  reported  as  the  EPTS.  However,  in  the  case 
of  RHA  (RC-37)  there  is  actually  a  gradual  change  in  the  stress  wave  at 
the  EPTS,  and  the  first  deviation  in  linearity  in  the  elastic  precursor 
could  be  an  indication  of  the  first  yielding  to  occur.  For  example,  the 
reported  EPTS  after  a  stress  wave  has  traveled  9. SI  mm  in  RHA  (RC-37)  is 
2.09  GPa,10  while  the  stress  corresponding  to  the  first  deviation  from  an 
approximately  linear  increase  in  stress  at  the  elastic  wave  front  is  1.77 
GPa.  A  proportional  adjustment  in  the  EPTS,  which  was  1.9  GPa  at  the 
region  of  maximum  fracture,  suggests  yielding  started  at  1.61  GPa  which 
is  less  than  the  nucleation  threshold  stress  of  1.65  GPa  -  a  reasonable 
result. 


V.  CONCLUSION 

A  nucleation  threshold  stress  cr  of  1.65  GPa  was  established  for 

no 

the  creation  of  cracks  in  rolled  homogeneous  steel  armor  (RC-37)  with 
impulsive  loading. 

The  failure  of  this  material  progresses  by  inclusion  cracking, 
followed  in  succession  by  the  growth  of  "hair-line"  cracks  away  from 
regions  of  broken  inclusions,  in-plane  crack  extension  from  inclusion  to 
inclusion  and  coalescence  of  non-planar  cracks  by  shear  on  connecting 
surfaces . 

The  nucleation  threshold  stress  of  1.65  GPa  corresponds  to  the 
initiation  of  the  "hair-line"  cracks. 

The  approach  used  to  establish  &no  depended  on  extrapolating  crack 

densities  developed  over  a  range  of  pressures  to  the  no-damage  condition. 
It  was  found  that  this  procedure  makes  sense,  but  the  tensile  stresses 
that  cause  failure  must  be  used  rather  than  the  impact  compressive 
stresses.  It  was  found  that  the  tensile  stresses  are  substantially  less 
than  the  compressive  stresses  because  of  elastic  wave  reflections  and 
unloading  at  void  interfaces.  These  factors  were  accounted  for  in  this 
work  with  computations  based  on  the  PUFF  computer  code  and  the  BFRACT 
subroutine. 
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